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Abstract 

Pulse shape analysis is an important background reduction and signal identification technique for next generation of 
76 Ge Ov/3/3 experiments. We present a study of the systematic uncertainties in one such parametric pulse-shape analysis 
technique for separating multi-site background from single-site signal events. We examined systematic uncertainties 
for events in full-energy gamma peaks (predominantly multi-site), double-escape peaks (predominantly single-site) 
and the Compton continuum near Qpp (which will be the dominant background for most Ov/3/3 searches). In short, we 
find total (statistical plus systematic) fractional uncertainties in the pulse shape cut survival probabilities of: 6.6%, 
1.5% and 3.8% for double-escape, continuum and y-ray events respectively. 
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1 1. Ov/3/3 in Germanium 

O 

The search for physics beyond the Standard Model has one of its most promising leads in neutrinoless double-beta 
decay (Ov/3/3) in particular, and neutrino physics more generally. Interest in Ov/3/3 is extremely well-motivated in the 
literature H] [2] [3] |4] [5] [5] . There are approximately 10 isotopes known to undergo two-neutrino double beta decay 
£C) (2v/3/3, i.e. 76 Ge — > 76 Se + 2<?~ + 2v), the slowest nuclear decay allowed in the Standard Model. These are also the 

isotopes of interest in the search for Ov/3/3 (i.e. 76 Ge — > 76 Se + 2e~ + Ov), a (3/3 mode forbidden by the Standard 
Model. If observed, Ov/3/3 would imply the existence of massive Majorana neutrinos[7] and could also lead to the 
discovery of other physics beyond the Standard Model. If T^J 2 could be measured in several isotopes, it would do 
{S) much to elucidate this physics |8 |. To make a meaningful comparison between T®J 2 measurements in several isotopes, 

the total experimental uncertainty must be sufficiently low. This places significant demands on the size of systematic 
uncertainties of the global experimental program. The importance of systematic uncertainties in potential results of 
Q\ Ov/3/3 searches motivates the investigation of systematic uncertainties in pulse shape analysis (PSA) presented here. 

76 Ge is one of the /3/3 isotopes under investigation in the current generation of Ov/3/3 searches (the most stringent 
T®J 2 limits already come from 76 Ge), and there are two next-generation 76 Ge experiments currently under develop- 
ment: Majorana||91 [TOl [TTl [T2l [131 and GERDA fBl . This article will present work performed in support of the 
Majorana experiment. Both Majorana and GERDA, as well as all Ov/3/3 searches (and indeed all searches for rare 
events), rely heavily on reducing backgrounds while retaining signals as effectively as possible. One of the ways 
that both experiments plan to reduce background and identify signal events is through PSA on the array of high-purity 
germanium (HPGe) detectors that will comprise the experiments. Because of the inherent complexity of making pulse 
shape cuts, systematic uncertainties in PSA-cut efficacy could be one of the leading contributors to the total systematic 
uncertainty budget of these experiments. An examination of these systematic uncertainties is the primary focus of this 
article. We will begin with a description of the detector used in this work, then move on to an overview of the use of 
PSA in HPGe detectors for Ov/3/3 searches up to and including this article. We will then detail our survey of systematic 
effects on PSA, and close with some discussion of the context of this work in the Majorana experiment. 
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2. Description of Experimental Apparatus 

The data presented in this article were taken using a CLOVER detector at Los Alamos National Laboratory. A 
CLOVER is a commercially available detector system from Canberra [15], consisting of four, n-type detectors in 
a single cryostat. The detectors in the CLOVER used in this study had two-fold azimuthal segmentation. These 
characteristics make the CLOVER a good off-the-shelf test bed for potential analysis techniques for Majorana. The 
individual detectors have a relative efficiency of 26%, this corresponding to a mass of roughly 750 g per detector. The 
CLOVER is instrumented with four high-resolution, cold-FET energy readouts (one for the central contact of each 
detector), and three low-resolution, warm-FET position readouts (corresponding to the left two, middle four and right 
two segments of the detector — see Figure [TJ. Coincidences between the energy and position readouts tell us which 




Figure I : A cartoon of the configuration and readout channels of the four detectors in the CLOVER detector used in this study. The central contact 
of each detector is numbered one through four, and left two, middle four, and right two segmentation contacts are labeled with an "L," "M," and 
"R" respectively. This Figure is adopted from Reference [ 15 ]. 

segment(s) recorded energy depositions. All four detectors share a common bias and preamplifier power. 

All contacts on the CLOVER were read out with a pair of DGF4C "Digital Gamma Finder" digitizer cards from X- 
Ray Instrument Associates (XIAI lfTBI . These boards operate under a Computer Automated Measurement and Control 
(CAMAC) architecture, and sample at 40 MHz (25 ns per sample) with a 14-bit analog to digital converter (ADC). 
Oversampling allows an effective resolution in the DGF4C boards of 16 bits. The data sets described in this paper fall 
into either "calibration sets" (used to calibrate the PSA algorithms) and "characterization sets" (used to characterize 
the efficacy of different event classes under PSA). For each individual systematic study, the characterization and 
calibration data sets were taken under identical hardware triggering and event selection conditions. 

3. Background Reduction Through Pulse Shape Analysis 

The lack of neutrinos in the final state of Qv/3/5 means that the only particles to carry away the available kinetic 
energy (Q^ = 2039.04 keV for 76 Ge fT7l ) are the two electrons. This pair of electrons will have a very limited 
range in germanium compared to backgrounds near Qpp, such as multiple Compton scatters from higher energy y 
rays. Most backgrounds near Qas w iU tend to produce interaction locations separated by up to several centimeters in 
germanium, allowing them to spread over a much larger volume and even across several detectors. events on the 
other hand will tend to stop in the HPGe detector in which they originated, depositing all their kinetic energy with 
high efficiency. Therefore, the signature of 0vf3j3 is a spatially localized (i.e. single-site) event with a well-defined 
energy (Q$?), uncorrected in time with any other event in the data stream. Pulse shape analysis, like most of the 
background tagging techniques employed in 76 Ge Ov/3/3 experiments, exploits the single-site nature of the signal and 
the multi-site nature of most backgrounds of similar energy. PSA tends to only be performed only on single-detector 
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events, because any multi-detector event is clearly multi-site. Nearly all PSA methods rely on digitized waveforms 
from the detector(s) (some techniques can be implemented only with analog electronics). The waveforms analyzed 
here are recorded from the output of the integrating preamplifiers of the detector system. We refer to such a waveform 
as a "charge pulse." Many PSA algorithms also analyze the "current pulse" of an event rather than the charge pulse. 
One can extract the current pulse by simply taking the time derivative of the charge pulse. 

Most PSA cut algorithms must first be calibrated on data containing known single and multi-site events, during 
which the calibration software associates typical waveform characteristics with each class of events (single and multi- 
site). For calibration data, we primarily used the full-energy y-ray peak at 1588 keV from 228 Ac for our population 
of multi-site events and the 1592 -keV double escape peak (DEP) of the 2614-keV y ray from 2()8 T1 for our population 
of single-site events. DEP events occur when a y-ray of sufficiently high energy undergoes electron-positron pair 
production in a detector, and both 51 1-keV y rays (created after the positron stops and annihilates) escape the detector 
without depositing any energy. DEP events are single-site because the only particles depositing energy in the detector 
of interest are charged leptons with a short range. Double escape peaks are described in more detail in ifTSl Chapter 
12]. Both the 1588-keV y-ray line and 1592-keV DEP are part of the 232 Th decay chain, making it operationally 
simple to collect both simultaneously. 232 Th has become something of a standard source in the examination of PSA in 
HPGe detectors because its decay chain has a relatively high-energy y ray at 2614 keV (giving a double-escape peak 
at 1592 keV), coupled with a spectrum otherwise free of strong lines above ~ 1700 keV. For some studies, we used 
the 1771-keV y ray and one of the several DEP peaks available in the spectrum of 56 Co. Although 56 Co has several 
high-energy y rays up to 3600 keV, giving DEP lines from 1600 to 2600 keV, spanning Q^g for 76 Ge, these y-ray 
lines also lead to a strong Compton continuum in this region, meaning that 56 Co DEP events can be hard to separate 
from continuum events for the purpose of PSA calibration. 56 Co also has a relatively short half-life (77.27 days), 
meaning that sources must be replaced rather frequently. For these reasons, 232 Th tends to be favored for calibration 
most PSA algorithms. Throughout this article we will refer to "DEP events" and "y-ray events" as those selected by 
making energy cuts on these features in either spectrum. The reverse of the techniques used in Ov/30 searches (i.e. 
rejecting single-site events while accepting multi-site events) could have potential applications to searches for rare 
nuclear decays to excited states, by allowing for the separation of weak y-ray lines from a strong continuum. 

3.1. Evolution of Parametric Pulse Shape Analysis 

PSA was used to varying degrees in both of the previous generation searches for Ov/3/3 in 76 Ge: the International 
Germanium Experiment (IGEX) lfl9ll20ll and the Heidelberg-Moscow experiment (HM) J2l] IZS |23] |24] |25) . PSA in 
the HM experiment is discussed in Reference ll25l Section 4.2]. The PSA work presented in this article is largely an 
extension of the PSA work done in support of the IGEX experiment, which is discussed at greater length in References 
[26, 27 , 28 1 . Our PSA uses parameters calculated from the charge and current pulses to quantify the single and multi- 
site nature of events from HPGe detectors. The way in which an algorithm associates single or multi-site events with 
numerical values of these parameters represents one of the key differences between parametric PSA algorithms. This 
is true of the changes in PSA from References Il26ll27l to [28 1 to that presented here. 

The IGEX experiment used three pulse shape parameters: 

• pulse width (defined as the time it took for the charge pulse to rise from 10% to 90% of its full amplitude), 

• front-back asymmetry (defined as the difference in the area of the first and second halves of the current pulse, 
normalized by the total area of the current pulse), and 

• normalized moment (defined essentially as the moment of inertia of the current pulse, were it treated as a mass 
distribution — see Reference 11281 for the exact form). 

The IGEX PSA algorithm was calibrated by creating a three-dimensional histogram of these parameters from a set of 
DEP events. The fraction of the total number of calibration events in each bin was then calculated. The bins of that 
histogram were rank ordered by that fraction, and stepping down the list, the cumulative fraction of the total number 
of events in the list was calculated for each element (so that the first element would have only the fraction of events 
in the most heavily-populated bin and the last element would have all of them). A cut was then made at an element 
on this list so that some nominal fraction (80% in most cases) of calibration events were preserved. This defined 
the signal space for the cut. When production data was later analyzed and its pulse shape parameters calculated, the 
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analysis software checked which bin in the three-dimensional PSA parameter histogram each pulse occupied. If that 
bin was in the signal space, it passed the cut. 

This method was modified for the preparation of Reference |28|. The detectors in the CLOVER are much smaller 
than those used in IGEX (< 0.75 kg versus > 2 kg). As a result, there as little difference in drift times from different 
regions of the individual CLOVER detectors, so the pulse width was largely degenerate between known single and 
multi-site events. Reference [281 therefore used only the second two parameters listed above. The calibration software 
created two-dimensional histograms for both y rays and DEPs from a 25%-subset of the characterization data. A fit 
was then performed to the parameter space histogram for each detector using a complicated analytical model (a sum 
of many two-dimensional Gaussian and Lorentzian functions). Both fit functions were then normalized so that their 
integral over the parameter space was one, allowing for their treatment as probability density functions (PDFs). When 
the full characterization data set was analyzed, the analysis software calculated the parameters for each pulse, and the 
PDF (fit function) for that detector was evaluated at those parameter values. If the PDF from the DEP (y ray) data had 
the higher value, then the event in question was classified as single-(multi-)site and passed (failed) the PSA cut. 

The PSA algorithm from Reference [28] was further modified for this work. Here, we calibrated the PSA algorithm 
with a data set completely independent of the characterization data (we did quantify the bias resulting from different 
amounts of overlap in between the calibration and characterization data — see Section 4.2 of this article for more 



detail). We also did away with the complicated fit described above. Instead, we simply normalized the PSA parameter 
histograms for DEP and y-ray events from each detector so that their integral was one, treating them as a PDF. The 
analysis of characterization data was then performed by checking the value of the PDF histograms corresponding to 
the pulse shape parameter values for each characterization event. 

Because a CLOVER detector is actually a four-detector array, we could select calibration DEPs in two different 
ways: one in which we made a simple energy cut in the single-detector spectra, and another in which we looked for 
two-detector events with the DEP energy in one detector, and one of the two e + annihilation y rays in another. The 
coincidence tag severely limits the number of DEPs available for calibration but also provides a very clean population. 
Most of the results presented in this article were made with single-detector, energy cut-selected DEPs. PSA calibration 
y-rays were always selected with a single detector cut. When both analyses are presented, the PDF constructed with 
non-coincidence tagged DEPs is denoted with a subscript "NC" (for "No Coincidence"). Figure [2] shows surface and 
contour plots of the PDF histograms for single-detector DEP and y-ray events. We could have also tagged DEPs by 
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Figure 2: A zoomed-in view of the PDF histograms from ' 
contour plots are also shown at top of both plots. 



2 Th data, populated with 1592-keV DEP (left) and 1588-keV y-ray (right) events. The 



demanding three-detector events: DEP in one detector, plus both e + annihilation y rays in other detectors. This would 
have provided an even cleaner population of calibration DEPs, but the efficiency for capturing such events in our array 
is impractically low for this study. 
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To evaluate the efficacy of different implementations of PSA in this article, we examine the ratio of the number 
of events which pass the PSA cut to the number of events on which the cut is performed, defining this as the survival 
probability. Survival probabilities were calculated for three different event classes: the DEP and y-ray events dis- 
cussed above and Compton continuum events near for 76 Ge. Understanding the survival probability of Compton 
continuum events is important because they represent the most likely backgrounds to experiments like Majorana. For 
continuum events, calculating the number of events in the pre and post-PSA cut spectra is simply an exercise in count- 
ing the number of events in an energy region. In this study, the region was either 2.0-2.08 MeV for the 232 Th runs 
(~ 40 keV on either side of Q3/3), or 2.04-2.08 MeV for the 56 Co source runs. The narrower energy region allowed 
the analysis to avoid the y-ray lines at 2015 and 2035 keV in the 56 Co spectrum. The statistical uncertainties are then 
just those arising from Poisson fluctuations. For peaks on top of continua, the situation is more complicated, because 
we need to separate out the strength of the peak from that of the continuum on which it rests. We do so by performing 
a standard^- 2 fit in the ROOT framework ll29l . In this case, the fit model is a flat background plus a Gaussian peak 
(or two, if multiple peaks reside near each other — this is the case for the 1588-keV y-ray peak and 1592 -keV DEP). 
To extract a peak's strength, we calculate its area from the fit parameters. The area uncertainty comes from that for 
the fit parameters and the expression for the area of a Gaussian using standard error propagation techniques. We can 
then take the ratio of these pre and post-PSA cut peak areas to calculate the survival probability for that cut. We again 
propagate uncertainty through the expression to obtain the survival probability uncertainty. Throughout this article, 
this is referred to as the "fit uncertainty." All uncertainties presented are one-sigma, rather than another confidence 
interval. 



3.2. Nominal Pulse Shape Analysis Efficacy and Comparison to Previous Work 

In Table [T] we present nominal values for the survival probability of pulse shape analysis cuts for y-ray, DEP, and 
continuum events along with total uncertainties (fit/statistical plus systematic). If the PSA algorithm worked perfectly 
and the different event classes had no overlap in the range of their pulse shape parameter distributions, the survival 
probability for DEP (y-ray) events would have a survival probability of 100% (0%). 



Table 1 : PSA Survival probabilities for different processes over a range of energy including systematic uncertainties detailed in Section[4] 



Process 


Energy (MeV) 


PSA Survival Probability (%) 


228 Ac y ray 


1.588 


29.0+1.9 


208 T j DE p 


1.592 


65.6 + 2.5 


Continuum 


2.0 - 2.08 


45.4 + 0.7 



For comparison, Reference 112811 found PSA performance for DEP, continuum and y-ray line events to be 75 ± 
2.0 + 2.1%, 43 + 0.9 + 3.0%, and 20 + 0.5 + 1.0%, respectively. The first uncertainty is statistical, and the second is 
systematic. Statistical uncertainties in Reference [28] were calculated assuming "root-N" Poisson fluctuations of the 
fit-peak area in each spectrum. This makes them somewhat smaller than the fit uncertainties presented in this article. In 
Reference [28], systematic uncertainties were also discussed in a rather narrowly-defined way. They came only from 
differences in the survival probabilities between data sets where the source was placed in different locations around 
the CLOVER cryostat. We estimated the magnitudes of uncertainties in PSA cut survival probabilities as a function 
of many systematic effects using the CLOVER. These were typically variations of the method used to implement the 
PSA cuts and are discussed at length in Section|4] 

It is clear that Reference [28] and the work presented in this article show a modest discrepancy in PSA perfor- 
mance. As stated above, the calibration data used for the PSA in Reference [28] was a 25%-subset of the character- 
ization data. This lead to a bias toward stronger separation of DEP and y-ray events. Furthermore, the complicated 
fit to the pulse shape parameter PDF in Reference |28 ] exacerbated the calibration overlap bias (which is discussed at 
greater length in Section 4.2 of this article). The large number of parameters in the PDF fit in Reference Il28l lead to 
artificial structure in the resulting fit functions. When we removed this overlap between the calibration and charac- 
terization data sets and re-ran the fit-based PSA during the preparation of this article, the survival probabilities come 
much closer to those presented in Table[T](24.3% for y-rays and 68.6% for DEP events). Note that this bias associated 
with the DEP and y-ray events is not present for the Compton continuum data, which was not used for calibration the 
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PSA in either algorithm. The continuum survival probability in this report and Reference [28 1 agree quite well, which 
is expected as the biases discussed above would only affect those data sets for which events within the calibration and 
characterization data sets are correlated. We do not present comparisons between this work and PSA from the IGEX 
or HM collaborations because the detectors used in those experiments were very different than those used to prepare 
this article and Reference 1 28 1 . 

4. Systematic Uncertainties in Pulse Shape Analysis 

We now present our detailed investigation of systematic uncertainties in the survival probabilities of different 
event classes under our PSA technique. This study can also likely serve as an explicit list of effects that should be 
investigated for any novel PSA technique. Many of these are effects which can be completely removed with the careful 
control of run conditions during the collection of calibration, characterization and production data (i.e. hardware 
triggering conditions). Some however, are true systematic uncertainties that cannot be removed operationally and can 
only be quantified (such as different calibration data sets or source location). 

4.1. PSA Parameter Histogram Resolution and Calibration Set Size 

The PDF histograms have a finite number of bins in each axis. We therefore examined how the binning of these 
axes affects both PSA performance and the required size of the calibration data set. To this end, we construct a series 
of plots showing PSA performance as a function of calibration set size for four PDF histogram axis binnings. Figure 
[3] shows this performance for: 10, 25, 50 and 100 bins per PSA parameter histogram axis. The range of the axes for 
each PDF histogram was held constant, so we are adjusting the resolution of the parameter space in each PDF. The 
calibration set size is characterized by the number of DEP events as they enter the data stream at a much lower rate 
than y-ray events. 

First, the figures show that as binning in the PSA parameter histograms becomes finer, the survival probabilities 
for each event class require larger calibration sets to reach their asymptotic values. Finer binning generally improves 
PSA performance, but there is a point of diminishing returns around 50 bins per axis. That is, increasing the number 
of bins from 10 to 25 to 50 improves separation of DEP and y-ray events, but increasing from 50 to 100 bins per axis 
does not, even with large calibration sets. Finer binning in PSA parameter space allows for the resolution of finer 
structure in the signal and background PDF histograms, however, once that structure is captured, increasing resolution 
does not seem to improve the algorithm. We therefore chose the optimal conditions for our ensuing PSA studies (i.e. 
the binning that provides the best performance with the lowest number of calibration events) of 50 bins per axis for all 
other studies presented. The asymptotically best DEP/y-ray separation is reached when there are > 500 DEP events 
per detector in the calibration data. This corresponds to between 2.8 and 5.6 million (11 million) total events with our 
experimental setup without (with) demanding the coincidence tag. 

4.2. Independence of Calibration and Characterization Data Sets 

As discussed above, it is important that the PSA calibration data be isolated from the data set used to characterize 
the efficacy of the cuts. If the data sets overlap, the cuts will appear more effective than they actually are because 
fluctuations in the pulse shape parameters for the different event classes will get frozen into the PDF histograms. This 
effect is most easily mitigated by completely separating calibration and characterization data. 

To quantify the bias in PSA performance due to overlapping calibration and characterization data sets, we used 
PDF histograms constructed from a single calibration set to analyze several characterization data sets of similar size 
(there were ten data files containing 7 x 10 5 events each). The overlap fraction is set by varying numbers of files in the 
characterization set in common with the calibration set and calculated by the ratio of the total size of the data files (in 
megabytes) in both the calibration and characterization data sets to the size of the characterization data set. We show 
the results of this study in Figure |4] 

The figure shows that as the overlap fraction increases, the separation of DEP and y-ray events increases for PSA 
trained on non-tagged DEPs. For comparison, we also included survival probabilities for PSA trained on coincidence- 
tagged DEPs. The term "overlap fraction" is something of a misnomer for the coincidence tagged DEPs because 
the calibration DEPs are all two-detector events and the characterization DEPs are all single-detector events (all 
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Figure 3: PSA efficacies as a function of the number of double-escape peak events per detector in the calibration data set for 10, 25, 50 and 100 
bins per PSA parameter axis. There are multiple data points for the small calibration set sizes because we constructed multiple calibration data 
sets. See the text for more detail. 
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Figure 4: PSA survival probabilities as a function of the fraction of the characterization data in common with the calibration data. 



calibration y-rays were single detector events regardless of how calibration DEPs were tagged), but we use the above- 
stated ratio of file sizes to calculate the overlap fraction for comparison. Figure|4]indicates that DEP and y-ray events 
for the non-tagged cuts are most strongly affected by sample bias, changing favorably 7.1% and 2.3%, respectively. 
As expected, the affect is essentially eliminated in the tagged cuts because the coincidence-tagged calibration DEPs 
never actually enter the characterization data stream. Finally, we see that the continuum survival probabilities are 
essentially independent of overlap fraction for the same reason: we do not actually train on continuum events, so there 
are none in the calibration data regardless of the overlap with the characterization data. 

4.3. Event rate and Pulse Pileup 

We now investigate the dependence of PSA performance on count rate in the CLOVER. Charge pulses from 
HPGe detectors have a short (typically < 1 /is) rise time, containing all of the information about energy deposition. 
In detectors with resistive feedback preamplifiers (such as the CLOVER), this is followed by a much longer fall time 
characterized by the preamplifier's RC time constant (typically ~ 50 //s). This fall time can be removed using a 
"Moving Window Deconvolution" |30|, but no such technique was employed in this work. Therefore, for our PSA to 
be maximally effective, the event rate had to be low enough that pulses from energy depositions in the detector are not 
deformed by the tail of the preceding pulse. If we assume an average time between pulses of ten to twenty preamplifier 
time constants, we estimate a maximum rate of ~ 1-2 kHz per detector. The pulse shape data presented in this article 
were taken over a range of rates from ~ 40 Hz up to ~ 900 Hz. Even at 2 kHz, (over twice our highest data rate), the 
fraction of the integral of a pulse overlapping another following by ten time constants (500 fis) should only be at the 
level of ~ 5 x 10~ 5 . For reference the single-detector rate expected in Majorana for production data is something like 
one count per hour, and no more than a few hundred Hz per detector for calibration. The principle source of dead time 
of our data acquisition (DAQ) system is actually communication across the backplane of the CAMAC crate. This is a 
much larger problem than pulse pile-up at rates approaching 1 kHz. The live time in the CAMAC DAQ system drops 
below 90% at around 80 - 90 Hz, drops below 80% by ~ 200 Hz and by 1 kHz has dropped to ~ 40%. Dead time 
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from backplane communication should have little or no effect on the discrimination power of PSA. Nevertheless, we 
looked for an effect. 

We adjusted the input count rate by placing different amounts of a 232 Th source two inches from the front face 
of the CLOVER within the lead shield. Event rates ranged from ~ 40 to ~ 500 Hz. At each rate, we collected 10.5 
million events for each of the three highest rates (98, 143 and 531 Hz), corresponding to ~ 1500 DEPs in the single- 
detector spectrum. The data set at the lowest rate (39 Hz) is twice as large corresponding to ~ 3000 DEP events. We 
therefore treat two halves of the low-rate data set independently. We trained the PSA cuts on the first two thirds of 
each data set and characterized their efficacy using the remainder. PSA results are reported in Table [2] 

Table 2: PSA Efficacies at different event rates. Survival probability uncertainties are just the error from the fit. Event rate uncertainties are the 
standard deviation of the rate reported by the DAQ software from run to run. 



Rate (Hz) 


DEP (%) 


y-ray (%) 


Continuum (%) 


530.8 + 39.1 


63.2 + 4.5 


31.1 ± 1.9 


46.8 ± 1.1 


143.1+0.1 


66.9 + 3.6 


29.2 ± 1.3 


46.4 ± 0.9 


97.6 + 0.3 


66.9 + 3.8 


33.1 ± 1.5 


48.2 ± 0.9 


38.8 + 0.1 


64.8 + 4.2 


35.9+1.7 


45.9 ± 0.9 


38.8 + 0.1 


68.1+4.0 


34.0 ± 1.7 


46.2 ± 0.9 



As expected, we find that the PSA survival probabilities depend quite weakly on the event rate. We do see that 
for the 531 Hz data, DEP survival drops slightly but not by a statistically significant amount. Still, because DAQ 
dead time from the CAMAC backplane is a problem at a considerably lower rate, we held the event rate in the rest of 
our studies down to below 100 - 150 Hz. Pulse shape distortion from pile-up could be a concern, however if a more 
modern DAQ system allowed for higher rates or in the case of two-step decays through states with ~ /is lifetimes. In 
this case, more sophisticated pulse processing would be required to remove the preamplifier's RC time constant. 



4.4. Different Calibration Sets 

Throughout the course of any experiment using PSA, experimenters will no doubt take many calibration data sets 
(assuming of course, that they use algorithms that require calibration data). Periodic calibration will also be necessary 
to check the stability of the detector array over time. We considered possible uncertainty in pulse shape analysis 
performance from calibration on different data sets taken under ostensibly the same conditions. 

The data sets described in Section 4.3 recorded at the three lowest rates (39, 98 and 143 Hz), were split into six 
equal parts, labeled A-F, providing 5 calibration data sets (A-E) and 1 characterization data set (F). There were a 
total of 60 700,000-event data files, so each of the six parts contained data files with roughly equal size taken at each 
event rate. We generated independent PDF histograms using each of the first five parts (A-E) and used those five PDF 
histograms to analyze the sixth data set. We compared the spread of the 5 deduced survival probabilities with their 
variance. If the spread is comparable to or smaller than the fit uncertainties for the individual results, then the fit errors 
place an upper limit on the calibration set uncertainty. If the spread is larger than the fit uncertainties, then we can use 
the difference between the spread and the fit uncertainty to estimate the additional systematic uncertainty arising from 
fluctuations in the calibration data. 

The last row of Table|3]shows us that the spread for DEP survival probabilities is markedly less than the average fit 
uncertainty. For y-ray and continuum events, it is only slightly larger, with an excess of 0.1% and 0.3% respectively. 
It is also worth noting that the spread in these values seem to be driven by calibration set D which has consistently 
higher survival probabilities for all event classes. We include these results as is, but it is probable that they are upper 
limits on this uncertainty. 



4.5. Triggering Conditions 

The DGF4C DAQ boards have three parameters that configure the triggering conditions: a threshold and two 
filtering times. According to Reference [ 16 1, one of the trigger times, called the "Flat Top Time" affects the digitizers' 
ability to trigger on slower or faster rising pulses and therefore should most strongly affect the population of single 
versus multi-site events. Multi-site events should tend to have slightly longer rise times than single-site events (though 
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Table 3: PSA performance on data set F, trained on data sets A through E. Uncertainties listed with the survival probabilities come from the peak 
fit. The row "Mean and SD" reports the mean and standard deviation of the previous five rows. The row "o"s o - cfjrft" is the difference between the 
standard deviation and the average fit uncertainty. See the text for more detail. 



Calibration 


Survival Probabilities from Analyzed Set 


Data Set 


DEP (%) 


y-ray (%) 


Continuum (%) 


A 


66.3 ± 2.7 


33.5 ± 1.1 


46.6 ± 0.6 


B 


65.8 ±2.7 


33.9 ± 1.1 


46.9 ± 0.6 


C 


66.0 ± 2.7 


33.8 ± 1.1 


47.3 ± 0.6 


D 


67.2 ±2.7 


35.9 ± 1.1 


48.2 ± 0.6 


E 


65.7 ±2.7 


32.6 ± 1.1 


45.9 ± 0.6 


Mean and SD 


66.2 ± 0.6 


33.9 ± 1.2 


47.0 ± 0.9 


(TSD-O" Fit 


-2.1 


+0.1 


+0.3 



the difference is quite small in the CLOVER). We therefore investigated the impact of this parameter by examining 
the data sets taken at 39 Hz for Section 4.3 (trigger flat top time = 0.2 yus) and collected another with the same 
configuration with a trigger flat top of 0.075 /-is. This spanned the range we had used for this filtering time in our 
entire experimental program with the CLOVER that still produces proper operation of the detector. We then analyzed 
it in the same way as in Section |4~3"1 (/.e. trained on two thirds of the data set and characterized the remaining third). 
We treated the two halves of the larger 0.2 fis data set independently. We compare the results in Table [4] To facilitate 
comparison, we also include the difference between the survival probabilities for data taken with each trigger flat top 
and each event class. 



Table 4: PSA cut survival probabilities comparing data taken under different triggering conditions. The 0.2-ps data was taken from the 39-Hz 
event rate data sets in Table [2] The survival probability uncertainties are once again those from the fit peak areas. The last two lines of the table 
(labeled "Difference") are the differences between the survival probabilities for each event class in the 0.075-/JS data set and each of the 0.2-/iS data 
sets. See the text for more detail. 



Flat Top 


DEP (%) 


y-ray (%) 


Continuum (%) 


0.075jus 


64.8 ± 3.7 


29.0 ± 1.5 


40.2 ± 0.7 


0.2//S 


64.8 ± 4.2 


35.9 ± 1.7 


45.9 ± 0.9 


0.2/js 


68.1 ±4.0 


34.0 ± 1.7 


46.2 ± 0.9 


Difference 


0.0 ±5.6 


-6.9 ±2.3 


-5.7 ± 1.1 




-3.3 ±5.4 


-5.0 ±2.3 


-6.0 ± 1.1 



Table|4]shows that the survival probabilities for DEP events remained unchanged when we lowered the trigger flat 
top time. Those for y-ray and continuum events however, are substantially lower with a shorter flat top time. This 
leads us conclude that PSA performance can be affected by triggering conditions. Trigger filter parameters should 
therefore be tuned for each detector in an experiment to maximize the efficacy of PSA as well as energy resolution 
and other characteristics. It is also important that pulse shape analyses of data sets with different triggering conditions 
be carefully considered. 

4.6. Analyzed Line Energy 

Figure [5] was constructed by calibration our PSA software on the 1576-keV DEP and 1771-keV y-ray lines from 
56 Co. We then plotted the survival probabilities for different lines in the spectrum as a function of energy using these 
calibration parameters. Note that the very low-energy y-ray lines in Figure[5] whose interaction is dominated by single 
photo-absorbtion (a very spatially localized process), have nearly identical behavior under PSA as the DEP lines, as 
expected. There is very little variation in the survival probability of y-ray lines above ~ 1 MeV Quantitatively, the PSA 
survival variation of DEP and y-ray events over 1 MeV has a standard deviation of 0.02% and 0.006% respectively. 
The average fit uncertainty is 0.06% for DEP and 0.005% for y-ray events. While the y-ray events show a variation 
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Figure 5: Line energy dependence of pulse shape analysis cuts for y and double-escape peak events. 



slightly in excess of the fit uncertainty, it is exceedingly small at the level of only 0.001%. Generally speaking, over 
the energy range available in 56 Co data (which more than covers the energies of interest for tagging Ov/?/? and 2v/3/3 
events), DEPs at different energies behave nearly the same, while y-ray events can interact through several different 
processes lfl8l . 

4.7. Calibration DEP Energy 

In principle, calibrating the PSA algorithm on one DEP line should be as good as doing so on any other (over the 
energy range accessible with a 56 Co source). We now investigate any variation in pulse shape analysis from the use 
of different DEP lines to train the cut. We generated the PDF histogram for each analysis with the 1771-keV y-ray 
line, and one of the five different DEP lines in the 56 Co spectrum from 1576 to 2429 keV. As with the preceding 56 Co 
study, DEP and y-ray survival probabilities are averaged for all twelve y-ray and five DEP lines above 1 MeV, and we 
quantify the excess uncertainty by looking at the difference between the standard deviation and average fit parameter 
uncertainty for the survival probabilities. We also include the survival probability for continuum events with energy 
from 2.04-2.08 MeV. We tabulate these results in Table [5] and find variation of the calibration DEP contributes an 
uncertainty of: 1.5% for y-ray events and no excess uncertainty for DEP or continuum events. 



4.8. Source Location Dependence 

Last, we examine the dependence of PSA performance on the position of the source used to collect the data. The 
source was placed above, in front of and to the side of the CLOVER. The source being placed above and to the side 
of the CLOVER changes the source's position with respect to the segmentation planes and the crystal axes of the 
detectors. Placing the source in front of the CLOVER illuminates the closed end cap parts of each detector as opposed 
to the parts of the detectors along their cylindrical axes. We return to generating PSA parameter PDFs using the 
1576-and 1771-keV DEP and y-ray lines from 56 Co. Survival probabilities for each event class are again averaged 
for all DEP and y-ray lines above 1 MeV and continuum events from 2.04-2.08 MeV, and we assess the uncertainty 
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Table 5: Average PSA survival probabilities for DEP and y-ray lines above 1 MeV and continuum events from 2.04-2.08 MeV for PSA trained on 
different DEP lines. <tsd is the standard deviation of the five survival probabilities in this table and <Xf„ is the average fit uncertainty for each cut. 



Calibration DEP 


DEP (%) 


r(%) 


Continuum (%) 


1576 keV 


62.2 ± 2.8 


29.4 ± 0.2 


50.0 ± 0.5 


2180keV 


62.0 + 2.8 


31.0 + 0.2 


49.9 ± 0.5 


2231 keV 


63.5 + 2.8 


30.8 ± 0.2 


50.4 ± 0.5 


2250 keV 


62.9 + 2.8 


31.8+0.2 


50.3 ±0.5 


2429 keV 


60.1 ±2.8 


34.1 ±0.2 


49.5 ± 0.5 


<J"SD - C Fit 


1.3-2.8 = -1.5 


1.7-0.2 = +1.5 


0.4-0.5 = -0.1 



from this variation again by looking at the difference between the standard deviation of the survival probabilities and 
uncertainty from the fit parameters. We tabulate these results in Table|6] Once again, we find no variation in the DEP 
or continuum survival in excess of the fit uncertainties, and a relatively small one for y-ray events at 0.8%. 



Table 6: Average PSA survival probabilities for DEP events and y-ray events above 1 MeV and continuum events from 2.04-2.08 MeV for runs 
with the source in different positions around the CLOVER, ctsd is the standard deviation of the survival probability and cr/r„ is the average fit 
uncertainty for each data set. 



Source Position 


DEP (%) 


y(%) 


Continuum (%) 


Above 


64.0 ± 6.8 


25.3 ±0.3 


45.7 + 0.8 


Front 


58.9 ±6.1 


27.3 ±0.3 


47.9 ± 0.9 


Side 


59.6 ± 7.6 


27.1 ±0.4 


48.2+1.2 


0~SD - V Fit 


2.8 - 6.8 = -4.0 


1.1-0.3 = +0.8 


0.2-1.0 = -0.8 



5. Discussion 

We now motivate the choice of specific values for the systematic uncertainty for the PSA probabilities. We 
summarize the results from Section|4]now in Table|7] Summing the effects in Table[7]in quadrature (thereby neglecting 
any potential correlations in these effects), we arrive at a total uncertainty of 1.9% for y rays and 2.5% for DEPs and 
0.7% for continuum events. Fractionally (that is, normalized to the survival probability for each event class), this 
comes to 6.6%, 3.8% and 1.5% respectively. Reference (8] sets a goal for total uncertainty budgets required to 
extract information about the mechanism underlying 0vj3/3 at roughly 20% to 50%, depending on the number of Ov/5y8 
measurements and model space of 0v/3/3 mechanisms considered. We can see that the fractional uncertainties in PSA 
performance will not be a major driver of that total uncertainty budget, but should nevertheless be kept as low as 
possible to increase confidence in the PSA algorithm. 

It is also important to note that the PSA performance demonstrated in this article used standard coaxial HPGe 
detectors read out with digitizers based on rather old technology. Current plans for data acquisition for the Majorana 
experiment include digitizers which sample at 100 MHz (as opposed to 40 MHz for the DGF4C boards used here). 
Majorana digitizers will also communicate over VME rather than CAMAC architecture, thus avoiding the backplane 
throughput problems encountered here. The Majorana collaboration also plans to field P-type Point Contact (PPC) 
detectors in its Demonstrator phase lfl3l . These detectors have been superior noise performance and excellent PSA 
discrimination power lFSTl [32l [331 . The Majorana collaboration has already begun purchasing "Broad-Energy Germa- 
nium" detectors (a particular style of PPC detector) from Canberra for deployment in the Demonstrator. The use of 
these detectors in lieu of standard semi -coaxial detectors as well as fielding more sophisticated digitization electronics 
should substantially improve PSA discrimination of the Majorana experiment over what has been shown here. 
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Table 7: Summary table of systematic effects impacting parametric pulse shape analysis with our CLOVER. The effects are quoted for PSA cuts 
trained on non-coincidence-tagged DEPs. Rows where the effect is marked with a "*" superscript can be held to zero operationally for production 
data taking. 



Effect 


Description 


y ray (%) 


DEP (%) 


Cont. (%) 


Calibration 


Different numbers of data files in common 


-2.6 x o-fi, 


+2.7 x o-fi, 


+0.5 x <T F i, 


overlap* 


between calibration and characterization data 








Event Rate* 


Rates < 143 Hz 











Different 


Trained on five independent data sets 


0.1 


— 


0.3 


Calib. Sets 


and analyzed same events 








Trigger 


Difference between two trigger settings 


3.5 




4.7 


Settings* 


minus fit uncertainty 








Analyzed 


DEP and y ray events above 1 MeV 


0.001 






Line 


from 56 Co data set 








Calibration DEP 


Trained on same y ray and different DEPs 


1.5 






Source Location 


Source above, in front, and to side of CLOVER 


0.8 








Systematic Sum 


1.7 




0.3 




Fit Parameter Uncertainty 


0.9 


2.5 


0.6 




Total Uncertainty 


1.9 


2.5 


0.7 




Nominal Performance 


29.0 ± 1.9 


65.6 + 2.5 


45.4 ± 0.7 




Fractional Uncertainty 


6.6 


3.8 


1.5 



6. Conclusions 

Much of this and previous work on PSA in germanium-based Ovfl/3 searches has focused on the notion of reducing 
backgrounds using PSA, but that is not the primary way in which the next generation of Ov/3/3 searches will exceed 
the sensitivity of the previous one. Enhanced sensitivity for future Qv/3/3 searches will be achieved primarily by 
building larger experiments out of cleaner materials, whereas novel analysis techniques will contribute to background 
reduction at a subdominant level. It is true however, that all analysis techniques employed in future Ovfi/3 searches 
will be examined for systematic uncertainties using a process similar to the one outlined in this article. This will 
be particularly important if several positive detections are made in several isotopes and the community undertakes a 
global Ov/3/5 analysis like the one in Reference (5). We used DEP and y-ray events as sample single and multi-site 
events, and quoting survival probabilities that are roughly a factor of two to three different between the two classes 
of events. It is far more likely that background events near Qpp will be Compton continuum events, and PSA (as well 
as most other background identification techniques) will not work as effectively on them. The T®J 2 sensitivity for an 
experiment is inversely proportional to the square root of the background level, given the same level of exposure(9j 
Section 4.6]. Reducing a background to 40% of its initial value (commensurate with the PSA efficacy for continuum 
events demonstrated here) would correspond to an increase in sensitivity of only around 60% ( Vl/0.4 = 1.6). Even 
if the combination of many techniques resulted in a factor of 10-20 reduction in background, this would correspond 
only to a roughly four-fold increase in T®J 2 sensitivity. While that level of impact on an experimental program is 
certainly not negligible, it is not the sole reason that the next generation 0v/3/3 searches should go to great lengths to 
tag backgrounds. Another extremely important reason to identify background events, is to demonstrate the single-site 
nature of any potential signal seen at or near Q^. This would do much to identify that peak as Ov/5j8 and not some 
hitherto unidentified background to a much higher degree of certainty. 
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